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ABSTRACT 

We study the problem of gravitational lensing by an isothermal elliptical 
density galaxy in the presence of a tidal perturbation. When the perturbation is 
fairly strong and oriented near the galaxy's minor axis, the lens can produce image 
configurations with six or even eight highly magnified images lying approximately 
on a circle. We classify the caustic structures in the model and identify the range 
of models that can produce such lenses. Sextuple and octuple lenses are likely 
to be rare because they require special lens configurations, but a full calculation 
of the likelihood will have to include both the existence of lenses with multiple 
lens galaxies and the strong magnification bias that affects sextuple and octuple 
lenses. At optical wavelengths these lenses would probably appear as partial or 
complete Einstein rings, but at radio wavelengths the individual images could 
probably be resolved. 



1. Introduction 

The first gravitational lens to be discovered, Q 0957+561, has a simple double image 
configuration (Walsh, Carswell & Weymann 1979). It was quickly followed by the first four- 
image lens, PG 1115+080 (Weymann et al. 1980). Together with Einstein ring lenses pro- 
duced by extended sources (e.g. MG 1131+0456, Hewitt et al. 1988), double and quadruple 
lenses nearly exhaust the list of configurations in the more than 50 known strong gravitational 
lenses.^ The only exceptions are 2016+112, which has three images and appears to require 
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the rare and complicated situation of two lens galaxies at different redshifts (Lawrence et 
al. f984; Nair & Garrett f997), and B f 933+503, which has ten images that can be nicely 
explained as a combination of three distinct sources, two of which are quadruply-imaged and 
one of which is doubly-imaged (Sykes et al. 1998). 

Theoretical studies have shown, though, that lenses with more than four images of a 
single source can exist. Schneider, Ehlers & Falco (1992) gave a mathematical analysis of 
caustic structures that yield additional images. Kochanek & Apostolakis (1988) surveyed 
models with two spherical galaxies at different redshifts and found that they can produce 
up to seven images. Witt & Mao (2000) found examples of models with an elliptical density 
galaxy and an external shear that can produce up to eight images, but they did not do a full 
survey of the models. Lenses with more than four images would be not only interesting to 
observe, but also very useful in determining the lensing mass distribution using the numerous 
constraints. 

In this paper we present a systematic classification of the caustics and image configura- 
tions for lenses consisting of an isothermal elliptical density galaxy and a tidal perturbation. 
We show examples of lensing cross sections and magnification distributions for different im- 
age numbers. In a forthcoming paper we will study the observability of lenses with double, 
quadruple, sextuple, and octuple image configurations. The outline of this paper is as follows. 
In §2 we review basic lens theory and the singular isothermal ellipsoid (SIE) and external 
shear lens models. In §3 we present an analytic classification of the caustics in a lens model 
with an SIE galaxy and a tidal perturbation approximated as an external shear. In §4 we 
study the stability of the caustics by adding a small core radius to the galaxy and letting 
the perturbation be produced by a neighboring galaxy or group. Finally, in §5 we discuss 
some of the observational consequences and applications of our results. 



The theory of gravitational lensing is discussed in detail by Schneider et al. (1992); here 
we summarize the features central to our analysis. The mapping between a source at angular 
position u on the sky and an image at angular position x is given by the lens equation, 



2. 



Methods 



2.1. Basic lens theory 



u — x — V(f)(x) , 




where the lensing potential is the projected gravitational potential of the lens mass. The 
potential is determined by the two-dimensional Poisson equation V 2 0(a ? ) = 2S(x)/S cr , where 
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£ is the projected surface mass distribution of the lens and the critical surface density for 
lensing is 

c 2 .Dos 

where D \ and D os are angular diameter distances from the observer to the lens and source, 
respectively, and D\ s is the angular diameter distance from the lens to the source. For a 
point source at u, there is an image at each root Xi of the lens equation ([[]). 



The brightnesses of the images are determined by the magnification tensor 

(3) 



*«- I, 



l-(f>, xx {x) -<P,xy{x) 

-4>,xy{x) 1 - 4>,yy{x) 



where subscripts denote partial differentiation, <p ti j = d 2 4>/dxidxj. The magnification of a 
point image at position x is given by | det Ai(x)\. In general there are several curves in the 
image plane along which the magnification tensor is singular and the magnification is infinite 
(det M.~ l = 0). These are called "critical curves," and they map to "caustics" in the source 
plane. Caustics mark discontinuities in the number of images, which leads to the key idea 
for our analysis: in order to determine the number of images produced by a lens model, it 
is sufficient to examine the caustics. If we know that a source arbitrarily far from the lens 
produces only one image, then we can imagine moving the source around and noting when it 
crosses a caustic to keep track of the total number of images. We discuss examples in §3.2. 



2.2. Model components 

We model a galaxy as a singular isothermal ellipsoid, because this model is not only 
analytically tractable but also consistent with models of individual lenses, lens statistics, 
stellar dynamics, and X-ray galaxies (e.g. Fabbiano 1989; Maoz & Rix 1993; Kochanek 1995, 
1996; Grogin & Narayan 1996; Rix et al. 1997). If we take the ellipsoid to be oblate with 
intrinsic axis ratio q% then its three-dimensional density distribution is 

= o 2 £3 1 U) 

9 27rGq 3 sin" 1 ^ R 2 + z 2 /q 2 ' 1 ' 



where G is the gravitational constant, a is the velocity dispersion and e 3 = J 1 — q 2 the 
eccentricity of the mass distribution, and (R, z) are usual cylindrical coordinates. For lensing 
we need the projected mass distribution. Choosing coordinates with the projected major 
axis along the x-axis, the projected surface density distribution is 

E 1 (5) 
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where bj = An — — (—} . (6) 

sin" 1 e 3 \cJ D os 1 ; 

The projected axis ratio q depends on the intrinsic axis ratio and the inclination angle i 
(where % = 0° is face-on and i = 90° is edge-on), 



q = \]ql sin i + cos 2 i . (7) 

The lensing properties of the isothermal ellipsoid have been given by Kassiola & Kovner 
(1993), Kormann, Schneider & Bartelmann (1994), and Keeton & Kochanek (1998). The 
lensing potential and the deflection angle (0 jX , 4> >y ) are 

(f> = x<i>,x + y<f>,y, (8) 

= V^W tan " {(^¥ X ) ' 
^ = 7^F tanh_1 (v5^7 y )- 

In the limit of a spherically symmetric mass distribution (a singular isothermal sphere, 
or SIS), the tangential critical curve is a circle with radius 

a\ 2 D, 



6sis = 4tt - -± (9) 

(in angular units); the maximum separation between images is ~ 2&sis- This "critical radius" 
therefore serves as a natural length scale in the lensing analysis, and we use it this way below. 
Typically 631s ~ 0.2" -3" for galaxy-scale lenses. 

In §3 we study models with a tidal perturbation approximated as an external shear. 
The potential and deflection angle for an external shear are 

= -i 7 r 2 cos2(fl-# 7 ), (10) 
(f) :X = — 7r cos(# — 2# 7 ) = —72; cos 2Q 1 — 73/ sin 29 1 , 
4>,y — 7 r sin(# — 29-y) = —7a; sin 2# 7 + jy cos 26* 7 , 

where 7 is the strength of the shear and # 7 is its direction angle, which is equal to the 
angle between the major axis of the galaxy and the shear axis because we have aligned the 
galaxy with the x-axis. Note that a shear described by (7, # 7 ) is equivalent to one described 
by (— 7, # 7 + 7r/2); we avoid this ambiguity by considering only shears with 7 > to be 
physical. 
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3. Galaxy+Shear Models 

A common lens system features a galaxy that is well modeled as an isothermal ellipsoid, 
plus a perturbation from objects at the same redshift as the lens galaxy (e.g. neighboring 
galaxies, or a group or cluster) or objects along the line of sight. For examples, see Hogg & 
Blandford (1994), Keeton, Kochanek & Seljak (1997), Kundic et al. (1997ab), Witt & Mao 
(1997), Tonry (1998), and Tonry & Kochanek (1999). To lowest order, the perturbation can 
be approximated as an external shear as in eq. fllPl). We begin by studying such galaxy+shear 
lens models because they admit a complete analytic treatment. In §4.2 we consider models 
in which the perturbation is instead produced by a discrete object like a galaxy or group. 



3.1. Critical curves and caustics 

To obtain a galaxy+shear lens model we merely add the lensing potentials in eqs. (§) 
and fllUD. The magnification from the joint model has a simple analytic form, 

det.M-^l-y'-v^ , l + 7Cos2(fl-0 7 ) 

r J{l + q 2)-( y l- q 2)cos29 



This lens model has a single critical curve, or curve of infinite magnification. From eq. ( O ) 
we can easily find a polar parametric form for this curve, 

r MW = f% l + 7COB2(0-* 7 ) 

1 ~T J(l + q 2 )-(l-q 2 )cos29 



The corresponding caustic is found by mapping the critical curve to the source plane with 
the lens equation ([I]). The caustic can be written in a cartesian parametric form, 

u cans (6) = [cos9 + 7 cos(# - 20 7 )] r crit (0) - 7y==f tan^ 1 (f cos#) , (13) 

^causW = [sin 6> - 7sin(6> - 20 7 )] r crit (0) &J tanh" 1 (f sin9) , 

V 1 - 9 



where ^ 



2(1 -g 2 ) 1/2 



l + g 2 )-(l-g 2 )cos2^ 



Note that in eq. ([I4) the parameter 6* is interpreted as the polar angle, but in eq. (0) 9 is 
simply a parameter that runs from to 2ir. 

The caustic is continuous but not smooth; in general it has four or more cusps. (See 
examples in §3.2.) We can find a simple equation that gives the location of the cusps. 
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Consider the parametric derivatives of the caustic, 



du cans /d6 
dv CSLVLS /d9 



-hi 



37sin2(#-# 7 ) + (1 + 7 cos 2(0 - # 7 ))£ 2 sin0cos0 



(i- 7 2)^r^ 



cos 9 + 7 cos( 



29„ 



sin 9 - 7 sin(# - 29 



7/ J 



(14) 



At a cusp the caustic stops and changes direction, so both du ca _ ns /d9 and dv caus /d9 vanish (and 
at least one of them changes sign). The only way for both derivatives to vanish simultaneously 
is for the common multiplicative factor in eq. ( |14| ) to vanish. We can simplify this to find 
that the condition for a cusp is 



7 sin 2 ( 



6, 



3(1 + q 2 ) -2(1 -q 2 ) cos 29 + (1 - q 2 )(sm 29 + 7 sin 2# 7 ) = . (15) 



If our model galaxy were non-singular, the lens model would have a second critical 
curve and caustic. Since it is singular, however, the second critical curve collapses to a point 
at the origin, and the corresponding caustic is considered to be a "pseudo-caustic." The 
pseudo-caustic has the cartesian parametric form 



"pscudo (0) 
L'pseudo (0) 



b 



tan 1 (£ cost 



(16) 



— x tanh 1 (£ sin 9) , 



where again < 9 < 2n, and £ is the same as in eq. (|i3|) . Note that the pseudo-caustic does 
not depend on the shear. 



3.2. Examples 

Figure 1 shows typical caustics and pseudo-caustics for galaxy+shear lens models, plot- 
ted using &sis from eq. @ as the natural length scale. In general the pseudo-caustic is 
smooth, while the true caustic has cusps that give it a diamond or astroid shape. As dis- 
cussed in §2.1, the caustics indicate where the number of images changes. A source outside 
both the caustic and the pseudo-caustic produces one image. When the source crosses to the 
inside of the pseudo-caustic it gains one additional image, which is faint, close to the galaxy 
center, and distorted radially relative to the galaxy; the pseudo-caustic is sometimes labeled 
the "radial caustic" because of the nature of the distortions. When the source crosses to the 
inside of the astroid caustic it gains two more images, which are bright, close to the corre- 
sponding critical curve, and distorted tangentially relative to the galaxy; hence this caustic 
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is sometimes labeled the "tangential caustic." The presence of the pseudo-caustic and the 
tangential caustic means that the number of images in these models is either one, two, or 
four depending on the location of the source. All of these features are generic to singular 
lens models that are not axisymmetric (see Schneider et al. 1992). 

The size and shape of the astroid caustic are determined by the net quadrupole moment 
of the lens model. When the galaxy and shear are aligned (# 7 = 0°, e.g. Figure la), the 
individual quadrupoles from the galaxy ellipticity and the shear combine to produce a larger 
astroid caustic than produced by either one alone. In other words, the area where sources 
produce 4-image lenses is larger. Mild misalignment between the ellipticity and shear twists 
the caustic (e.g. Figure lb). As the shear becomes more misaligned (e.g. Figure lc), the 
quadupole from the shear partially cancels the quadrupole from the ellipticity; this is why 
Figure lc has the smallest astroid caustic despite having the largest shear. We return to this 
point in the discussion (§5). 

Figure 2 shows that when the shear is nearly orthogonal to the ellipticity (# 7 ~ 90°), 
it can nearly cancel the effects of the ellipticity to produce a caustic that is quite small. 
Nevertheless, the caustic is very interesting because it shows qualitatively new features not 
seen in Figure 1. (Similar examples can be found in Witt & Mao 2000.) The caustic folds over 
on itself in features called "swallowtails" (see Schneider et al. 1992). Because the number 
of images increases by two when the source crosses the caustic, a source inside a swallowtail 
produces six images (e.g. Figure 2b). The swallowtails are sensitive to the angle between the 
ellipticity and shear. They grow larger as the misalignment increases, until with near perfect 
misalignment the swallowtails overlap (e.g. Figure 2c). The region of the source plane where 
swallowtails overlap corresponds to sources that produce eight images. In other words, the 
interaction of the ellipticity and shear makes it possible to have new image configurations 
with six or even eight images. 

Not all combinations of ellipticity and shear can produce these new image configurations. 
Figure 3 shows caustics for three cases where the shear is orthogonal to the ellipticity (# 7 = 
90°). If the shear is small (Figure 3a), the system is dominated by the ellipticity and any 
swallowtails that exist are small. If the shear is large (Figure 3c), the system is dominated by 
the shear and again any swallowtails are small. For a given ellipticity, only a narrow range 
of orthogonal shears can produce overlapping swallowtails (Figure 3b). In §3.3 we discuss in 
detail the range of models that can produce 6- and 8-image configurations. 



- 8- 



3.3. Models that can produce more than 4 images 

Appendices A and B give a mathematical analysis of the range of models that have 
swallowtails and thus can produce 6- or 8-image lenses; we summarize the results here. We 
use the presence of swallowtails in the caustic to indicate that a model can produce more 
than four images. This approach does not directly indicate the probability of observing a lens 
with more than four images. Estimating this probability requires detailed computations of 
lensing cross sections and magnification distributions for realistic sets of lens environments. 
We discuss these issues briefly in §5 and plan to study them in more detail in a forthcoming 
paper. For now we seek to delineate the conditions under which a lens can produce more 
than four images. 

Figures 4-6 show the envelope of swallowtail models in the (q, 7) plane for different 
values of the shear angle 6 1 . When the shear is orthogonal to the ellipticity (Figure 4), 
the envelope encloses small shears for modestly flattened galaxies (q <; 1) and larger shears 
for more flattened galaxies (q <C 1). Swallowtails can exist for 7 arbitrarily small and q 
arbitrarily close to unity (ellipticity arbitrarily small) , provided that the shear and ellipticity 
are in the narrow band where they properly balance each other. A relatively broad band 
of the swallowtail models also have overlapping swallowtails and thus can produce 8-image 
lenses. 

If # 7 changes by even a few degrees away from perfect misalignment, however, the enve- 
lope pinches off at the high-g/low-7 end (Figure 5). In other words, obtaining swallowtails 
requires more ellipticity and shear. Also, the range of models with overlapping swallowtails 
quickly disappears (not shown). From these results we conclude that 6- and 8-image lenses 
produced by small tidal perturbations (e.g. 7 ~ 0.1) are likely to be very rare because they 
require a special combination of galaxy axis ratio and shear misalignment angle. However, 
with stronger perturbations the range of swallowtail models is considerably larger. Thus we 
predict that even though 6- or 8-image lenses may still be rare, they are more likely to occur 
when the perturbation is strong (such as a second galaxy near the primary lens galaxy) than 
when the perturbation is weak. Since strong perturbations may not be well approximated by 
the external shear model, in §4.2 we study models using a full treatment of a perturbation 
from a nearby galaxy or group. 

The mathematical analysis does reveal that swallowtails can exist for smaller shear 
misalignments; while the results are formally interesting, they are physically implausible 
because they require highly flattened galaxies and very large shears (Figure 6). For 9 y ~ 
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60°, there are two envelopes that are comparable in size.0 For 45° < 9 1 < 50.45° (0 > 
cos2# 7 > — 1/\/28, see eq. p^§| ), one of the envelopes stretches to arbitrarily large shears. 
For 39.55° < 8 y < 45° (1/V28 > cos 26> 7 > 0), the finite envelope disappears but the infinite 
one remains. This envelope finally disappears for 9 7 < 39.55° (cos26* 7 > l/\/28). With the 
strong perturbations required by these envelopes, however, the simple shear approximation 
almost certainly breaks down. Thus this analysis is probably not fully valid, but it does 
suggest features to look for when studying models with very strong perturbations, such as 
interacting galaxies. 

4. Stability of the Caustics 

4.1. Non-singular lens models 

Although we have studied singular lens models for analytic convenience, real galaxies 
may have a small but finite core radius. It is important to understand whether swallowtails 
are robust under the addition of a core radius. Adding a core radius adds one faint image 
near the center of the galaxy for any configuration with more than one image, so the total 
number of images is always odd (see Schneider et al. 1992). Figure 7 shows generalizations 
of the caustics in Figures 2b and 2c to a finite core radius s, using the lensing properties of a 
softened isothermal ellipsoid given by Keeton & Kochanek (1998). The swallowtails shrink 
as s increases - probably because when we increase s with fesis held fixed, we decrease the 
galaxy mass and hence reduce the contribution of the ellipticity to the potential. However, 
no convincing case of a faint central image has been observed, and this limits the core 
radius (s/&sis ^ 0.1, see Wallington & Narayan 1993; Kochanek 1996) and suggests that 
lens galaxies are quite cuspy. The presence of such small cores would not significantly affect 
our results. 



4.2. Two-galaxy models 

As discussed in §3.3, the galaxy+shear models suggest that the perturbation must be 
relatively strong in order to produce swallowtails. For such strong perturbations, the external 
shear approximation may not be justified, so in this section we examine a simple model in 
which the perturbation is produced by a second mass distribution representing a neighboring 



2 The second envelope exists even for larger misalignments, as indicated by Figure 12 in Appendix A. 
However, it is too small to be seen in Figure 5. 



-10- 



galaxy or group. We again use a singular isothermal mass distribution for the perturber, 
but for simplicity we assume it is spherical. 

The perturber is described by three physical parameters: its velocity dispersion 02 (or 
equivalently its critical radius &sis,2 given by eq. |9|),f] and its position relative to the lens 
galaxy, given by the projected distance d from the galaxy center and the angle Q 1 from 
the lens galaxy's major axis. To generalize the galaxy+shear models, it is convenient to 
characterize the perturber not by its velocity dispersion o<i but rather by the strength of 
the perturbation. To lowest order, the perturbation is equivalent to an external shear with 
strength 

W (2d) 

7 1 - b sis , 2 /(2d) ■ 1 ) 

It is important to understand what this strength means. All perturbers with a given strength 
7 are equivalent to each other and to an external shear of strength 7 - to second order in 
the potential. The differences enter only in terms of third order and higher. Nevertheless, 
we show here that the differences are important. 

We note that since the perturbation strength 7 is a combination of the velocity dispersion 
and distance, perturbers with a given strength but different distances must have different 
velocity dispersions, 

'0-2 \ 2 &SIS.2 d 7 



- 



(1* 



vCi/ &sis,i &sis,i 1 + 7 
In other words, the more distant the perturber is, the more massive it has to be, and vice 
versa. 

The top panels in Figure 8 show the caustics for a galaxy+shear model with q = 0.5 
and # 7 = 90°. They are similar to Figure 3 but include more values of 7. The other panels 
in Figure 8 show the generalization to two-galaxy models, using the same values of 7 and 
distances d/bsis,i = 5,10,20,40. Note that the galaxy+shear models in the top panel are 
equivalent to two-galaxy models with d — > 00. Figure 9 is similar, but shows the caustic 
structures for 6 1 = 88°. From eq. (|18"D, the models in Figures 8 and 9 have perturbers that 
range in mass from the scale of a galaxy to that of a cluster. 

The examples show two important effects. First, the swallowtails tend to be larger in 
two-galaxy models than in equivalent galaxy+shear models, and they grow as the distance d 
decreases. The various models differ only in higher order terms, but those terms apparently 



3 In this section we use a subscript 1 to denote the main lens galaxy and a subscript 2 to denote the 
perturber. We continue to use the critical radius &sis,i of the main lens galaxy as the natural scale length. 
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strengthen the perturber's effects. Equivalently, the range of models that produce swal- 
lowtails is larger. Second, in galaxy+shear models the symmetry of the shear implies that 
there are always two identical swallowtails. By contrast, in two-galaxy models one of the 
swallowtails is often quite large, while the other is either small or absent. As a result, the 
area in the source plane where swallowtails overlap is small or absent, so the models have a 
limited ability to produce 8-image lenses. 

We conclude from these examples that the qualitative features we saw in the galaxy+shear 
models are stable when we change the source of the perturbation. In fact, the range of swal- 
lowtail models seems to increase with more realistic treatments of strong or close pertur- 
bations. However, it appears that the overlapping swallowtails required to produce octuple 
lenses are rare and not very stable. 

5. Summary and Discussion 

We have studied the lensing properties of an isothermal elliptical density galaxy in the 
presence of a tidal perturbation to classify the caustic structures and identify different image 
configurations. In most cases the models have a simple caustic structure corresponding to 
standard 2-image and 4-image lens configurations. However, when the tidal shear has a 
magnitude and direction appropriate to partially cancel the effects of the galaxy's ellipticity, 
the caustics develop complicated swallowtail features that correspond to 6- and 8-image 
configurations that have not yet been observed. We gave a complete analytic treatment 
of the case of a singular galaxy with a perturbation modeled as an external shear, but we 
showed that the caustic structures are stable when one adds a small core radius or uses a 
more realistic treatment of the perturbation. In fact, the swallowtail caustic structures are 
generally bigger with the more realistic perturbation than with the shear approximation, 
which may enhance the likelihood of observing lenses with more than four images. 

Our analysis has several observational consequences. First, sextuple and octuple lenses 
are likely to be rare because they require special lens configurations. In fact, they will 
probably be found only when there is a second galaxy close enough to the lens galaxy to 
provide a strong perturbation. While this situation is uncommon, it is not exceedingly rare: 
at least five lensesQ appear to have multiple lens galaxies within the Einstein radius. Still, 
to quantify the likelihood of finding a sextuple or octuple lens, it is necessary to compute 



4 B 1127+385 (Koopmans et al. 1999), B 1359+154 (Rusin et al. 1999), B 1608+656 (Koopmans & 
Fassnacht 1999), 2016+112 (see Nair & Garrett 1997 and references therein), and B 2114+022 (Augusto, 
Wilkinson & Browne 1996; Jackson et al. 1998). 
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the lensing cross sections for various combinations of the galaxy axis ratio q and the shear 
amplitude 7 and misalignment angle # 7 . Figure 10 shows cross sections for q = 0.5, 7 = 0.22, 
and three values of # 7 that correspond to the cases shown in Figure 2. The configurations with 
more images have smaller cross sections and thus smaller probabilities for being observed. 
However they also have significantly higher magnifications, which means that magnification 
bias will be important (e.g. Turner 1980; Turner, Ostriker & Gott 1984). Magnification 
bias will mitigate the effects of small cross sections to increase the likelihood of observing a 
sextuple or octuple lens. Clearly a realistic prediction of the lensing probabilities will require 
computation of cross sections and magnification bias for many combinations of (q, 7, 6> 7 ) 
weighted by realistic populations of lens galaxies and perturbers. We plan to address these 
issues of observability in a forthcoming paper, both in the context of current surveys and 
of the Next Generation Space Telescope, where large numbers of lenses are expected (e.g. 
Barkana k Loeb 1999). 

Second, if discovered the sextuple and octuple lenses that we have described will be 
easy to identify because the lensed images lie approximately on a circle (see Figure 2). Any 
lens with more than four images that do not trace a circle must not be of the type described 
here. Indeed, in B 1933+503 there are ten images not on a circle, and it is thought that 
they must be associated with three different sources (Sykes et al. 1998). In B 1359+154 
there are six radio sources with four sources in a standard quadruple lens configuration, plus 
two sources inside the configuration whose interpretation is not clear (Myers et al. 1999). 
Although observations and models suggest that there may be multiple lens galaxies (Rusin 
et al. 1999), the fact that the six sources do not follow a circle suggests that lens cannot be 
explained by swallowtails produced by the lens galaxies. 

Third, observed sextuple or octuple lenses would be very useful for constraining models 
of the lensing mass distribution. The "extra" images beyond the standard two or four would 
provide additional position and flux constraints. Even better, the sensitivity of the caustic 
structures to the lens galaxy ellipticity and the shear amplitude and misalignment angle 
means that the mere existence of six or more images should place strong constraints on 
those properties of the model. As a result, a sextuple or octuple lens could represent a 
wonderful ability to break the common degeneracy between the lens galaxy shape and the 
shear from the surrounding environment (see Keeton et al. 2000a). Such a lens might even 
serve as the long-sought "golden lens" for measuring the Hubble constant H (see Schechter 
2000) provided that a time delay could be measured, although the time delay might be 
relatively short (a few days) because the images are all close to the critical curve. 

These applications are predicated on the ability to resolve the individual images, and 
this ability might be limited when we include the effects of finite source size. Figure 11 
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shows that the source does not have to be very large before the images smear into a partial 
or complete Einstein ring. This happens for -R src /&sis ~ 0.01, where typically &sis ~ 1"— 3" 
for galaxy-scale lenses. Radio surveys can achieve sub-milli-arcsecond resolution with VLBI 
or VLB A mapping (e.g. the CLASS survey, see Browne 2000 and references therein), so 
they should still be able to resolve individual images and thus find lenses amenable to these 
applications. By contrast, optical surveys (such as the Sloan Digital Sky Survey, see Gunn 
et al. 1998 and Fischer et al. 1999, or the Next Generation Space Telescope, see Barkana 
& Loeb 1999) are more likely find partial or complete Einstein rings. Still, new techniques 
for modeling Einstein rings show that they are very useful for constraining not only the lens 
model but also the intrinsic shape of the source (Keeton, Kochanek & McLeod 2000b). 

Finally, our study may be relevant for the so-called ellipticity "crisis," where lens galax- 
ies are inferred to have larger ellipticities than the observed early-type galaxies (see Kochanek 
1996 and references therein). From Figure 1, it is clear that the caustic structures are the 
largest when the shear is aligned with the major axis of the lensing galaxy so the ellipticity 
and shear act coherently; conversely, the caustic structures are the smallest when the shear 
is orthogonal to the major axis since the shear partially cancels the lens ellipticity. Obser- 
vationally, this means that for a sample of quadruple lenses the shear may preferentially 
lie along the galaxy major axis; a good example is B 1422+231 (Hogg & Blandford 1994). 
The presence of an aligned shear is difficult to infer from lens models due to a degeneracy 
in the lens equation (see Witt 1996): naive lens models simply imply a model ellipticity 
equivalent to the combination of the true ellipticity and the shear. This effect may generate 
a bias toward larger inferred ellipticities for lens galaxies, although a full account of this bias 
awaits further investigations. 

Acknowledgements. We thank Peter Schneider for helpful discussions, and the anony- 
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A. Models that can produce at least 6 images 

In this Appendix we compute the range of galaxy+shear models (see §3) that can 
produce at least six images. Because the caustics determine the image number, finding 
models that can yield at least six images is equivalent to finding models in which the caustic 
has swallowtails. We saw in §3.2 that in models without swallowtails the caustic has four 
cusps, while in models with swallowtails the caustic has more than four cusps. Thus to 
identify models with swallowtails it is sufficient to find models that have more than the 



standard four solutions to the cusp equation (|T5|). 



-14- 



In other words, the envelope bounding the region in parameter space where models have 
swallowtails is located where the cusp equation develops additional pairs of solutions. If the 
cusp equation (|I5D is written as f{9) = 0, the place where additional solutions appear is 
defined by 

df(9) 



f(0) = 



and 



39 



0. 



(Al) 



These equations yield two polynomials in sin26 l and cos 2$, so we can use the resultant 
method (e.g. Walker 1955; Erdl & Schneider 1993) to eliminate 9. We find that the envelope 
is given by roots of the equation 



F(7,0 7 ,g) = ^a;(# 7 ,g)(l-g : 

i=0 



2^6-^ i 







where the coefficients a 



q) are: 



(A2) 



a 
a i 

«2 

a i 

a 5 
o 6 



18c 7 (l +g" 



(1 + g 4 )(74 + 49c 2 ) - 40g 2 + 334c 2 g 2 

12c 7 (l + q 2 ) [(1 + g 4 )(29 + 4c 2 ) + 50g 2 + 64cV 



16c T (l - g 2 ) 2 (l + 34g 2 + q A ) + 2c 2 (7 + 58g 2 + 7g 4 )(43 + 22g 2 + 43g 4 ) 



- 7 \ i / \ i ""3 ' i j • - 7 

-3(1 + 380g 2 - 2058g 4 + 380g 6 + q 



6c 7 (l + q 2 ) [4c 2 (l - g 2 ) 2 (19 + 358g 2 + 19g 4 ) - 1 - 500g 2 + 4890g 4 - 500g b - g b 
128c 4 (l - g 2 ) 4 (l + 34g 2 + q 4 ) - 3c 2 (l - g 2 ) 2 (l + 644g 2 - 6906g 4 + 644g 6 + g 8 ) 
+36g 2 (l + 34g 2 + g 4 ) 2 (A3) 



with c 7 = cos26' 7 . For any set of parameters ( n f,9 1 ,q), F > means that the caustic does 
not have swallowtails, while F < means that the caustic does have swallowtails. Thus each 
pair of solutions to F = gives an envelope bounding a region in parameter space in which 
swallowtails occur. 

For 6L = 0° and 90° we have c 7 = +1 and —1 and the envelope equation simplifies to 



F(rf,9^q)= l-q 2 ± 1 (5 + q 2 ) l-q 2 ± 1 (l + 5q 



3 r 



(A4) 



With c 7 = +1 there are no solutions, and hence no models with swallowtails. With c 7 = — 1 
the envelope of swallowtail models is 



1-g 2 1-g 2 

< 7 < 



5 + q 2 



l + 5q 2 



(A5) 
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More generally, for given values of # 7 and q the envelope equation (A2) is a 6th order 
polynomial in 7 whose roots are easy to find numerically. Before using a root finder, however, 
we should understand how many roots to expect and what their general ranges are. This 
requires that we examine how the polynomial F depends on location in the (q, # 7 ) plane. 
First consider how the number of roots changes as q varies. Additional pairs of roots appear 
when F = and dF/dq = 0, which can be combined using the resultant method to eliminate 
7 and obtain the condition 



G(cy,q) 



32(1 -g 2 ) 2 (l 



970g 2 + g 4 )c 4 - (7 - 986g 2 + 7q 



4\2 



+ 



17(1 + q 8 ) + 60196(g 2 + q 6 ) + 824358g 4 



where again c 7 = cos 26 1 . This is a second order polynomial in c 2 so its solution is 



(A6) 



1 



64(l-g 2 ) 2 (l + 970g 2 + g 4 ) 



17(1 + q 8 ) + 60196(g 2 + q 6 ) + 824358g 4 



+81(l + g 2 Wl + 322g 2 + g 4 



(A7) 



(The second solution of the quadratic equation is unphysical because it has c 2 < 0.) The 
right-hand side of eq. (|A7|) is between and 1 for all < q < 1, so there is always a physical 
solution; in fact, there are two solutions, one for c 7 > and one for c 7 < 0. The two solutions 



intersect when c 7 = 0, which occurs at q — y 7/(493 + 90y30) = 0.0843. Next, consider the 
behavior of the polynomial F for large 7, which is controlled by the coefficient ag. If > 
then F > for large 7, so there are no swallowtails for large 7. However, if < then 
F < for large 7, so there are swallowtails for arbitrarily large 7. Thus the behavior of the 
polynomial changes qualitatively when a 6 changes sign. Since a 6 is a quadratic polynomial 
in c 2 (see eq. |A3|) , it is easy to find that a 6 = if 



256(1 -g 2 ) 2 (l + 34g 2 + g 4 ) 



[1 + 6Uq 2 - 6906g 4 + 644<f + q 8 



±(1 - 253<f - 253g 4 + q b ) Jl - 254g 2 + q 



(A8) 



The right-hand side is real for |g| < 8 — 3y/7 = 0.0627, and at this point it has c 7 = ±l/\/28 
or 9 7 = 39.55° and 50.45°. 

With these results we can understand the (q, 7 ) plane as shown in Figure 12. The curves 
given by eqs. (|A7|) and (|A^) define nine different regions where the swallowtail envelopes (the 
parameter regions with F < 0) have the following properties: 



none : no envelopes 



- 16 - 





1+ 


one envelope, with 7 > 




2+ 


two envelopes, both with 7 > 


1+, 


1- 


two envelopes, one with 7 > and one with 7 < 


1+ 


, oo 


one finite envelope with 7 > 0, and envelopes with 7 




1- 


one envelope, with 7 < 




2- 


two envelopes, both with 7 < 


1-, 


1+ 


two envelopes, one with 7 < and one with 7 > 


1- 


, oo 


one finite envelope with 7 < 0, and envelopes with 7 



Negative shear is unphysical (see §2.2), so we are interested only in envelopes with 7 > 0, 
which occur in six of the nine regions: (1+); (2+); (1+, 1— ); (1+, 00); (1 — , 1+); and (1 — , 00). 

With this detailed knowledge of the (q, 7 ) plane, we can use a numerical root finder with 
eq. ( |A2| ) to obtain the envelope of swallowtail models. The results are shown and discussed 
in §3.3. In particular, Figures 4-6 show envelopes in the (9,7) plane for different values of 
6 1 . Each envelope corresponds to a particular horizontal line in Figure 11. Each point where 
this line crosses a curve in the (q, # 7 ) plane corresponds to a cusp in the envelope in the 
(<?, 7) Plane. 



We can also ask what range of swallowtail models produce overlapping swallowtails 
that bound 8-image regions (e.g. Figure 3b). We cannot answer this question analytically 
for arbitrary values of the shear angle # 7 . However, from the examples in §3.2 we expect that 
overlapping swallowtails occur only when the shear is nearly orthogonal to the galaxy, and 
the case with 7 = 90° can be studied analytically. When 6* 7 = 90° the system has reflection 
symmetry, so the points on the caustic with 6 = and 6 = n/2 are always cusps, no matter 
how convoluted the rest of the curve is. (These cusps are indicated in Figure 3.) Moreover, 
the 6 = cusp always opens to the left, and the 6 = n/2 cusp always opens downward. Label 
the 6 = cusp position (uh,0) and the 6 = tt/2 cusp position (0, vy) (H for horizontal, V 
for vertical). From Figure 3 we see that there are overlapping swallowtails only if uh and 
vy are both positive. (If one is positive and one negative then there are swallowtails that do 
not overlap.) We can use eq. (|13|) to rewrite the conditions uh > as follows: 



B. Models that can produce 8 images 




(Bl) 
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We can similarly rewrite the condition vy > 0: 



v v >0 <=► 7 >7v(g) = f44-T' 



where E v (q) = tanh 1 — ^ 

V 1 — V v 

For any galaxy axis ratio < q < 1, these functions satisfy 7y(<?) < 7#(<?)- Thus the 
condition for overlapping swallowtails is 7y(g) < 7 < 7h(<?), and this result is shown in 
Figure 4. 
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Fig. 1. — Sample caustics and pseudo-caustics for models with a singular isothermal ellipsoid 
galaxy and an external shear. The smooth elongated curves are the pseudo-caustics and the 
diamond-shaped curves are the astroid caustics. All panels have a galaxy with a projected 
axis ratio q = 0.5 and a shear whose magnitude 7 and direction 6 1 are indicated. The axes 
are labeled in terms of the natural lensing length scale 631s from eq. 
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Fig. 2. — Examples of "swallowtail" caustics produced when the ellipticity and shear are 
nearly orthogonal. The left-hand panels show the tangential caustics and sample source 
positions; the pseudo-caustics are larger than the frames. The right-hand panels show the 
corresponding tangential critical curves and image positions. All models have q = 0.5 and 
7 = 0.22 and the specified shear direction 9 1 . The points show that a source inside the 
astroid but outside the swallowtails produces 4 images (case a); a source inside a swallowtail 
produces 6 images (case b); and a source inside overlapping swallowtails produces 8 images 
(case c). The total magnification of the sample images in cases a, b, and c is 76.8, 178.3, 
and 266.2, respectively. 
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Fig. 3. — The dependence of swallowtails on the shear. All models have q = 0.5 and an 
orthogonal shear {6 1 = 90°) with the specified magnitude 7. The points indicate the cusps 
with 9 = (on the horizontal axis) and 6 = n/2 (on the vertical axis); see Appendix B for 
details. 
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0.2 0.4 0.6 0.8 1 

axis ratio q 

Fig. 4. — The envelope of galaxy+shear models that can produce at least 6 images in the 
(q, 7) plane, for a shear orthogonal to the lens galaxy (0 7 = 90°). Models in the cross-hatched 
region have overlapping swallowtails and can produce up to 8 images. Models in the shaded 
region have non-overlapping swallowtails and can produce up to 6 images. Models outside 
the shaded region do not have swallowtails and can produce at most 4 images. The curves 
bounding these regions are given by eqs. (|A5|) , flBl|), and ( p2|) in Appendices A and B. The 
three filled points indicate the locations of the three sample models shown in Figure 3. 



-24- 




0.2 0.4 0.6 0.8 1 



axis ratio q 

Fig. 5. — Similar to Figure 4, but for different values of the shear angle near 9 1 ~ 90°. Only 
the swallowtail envelope is shown (not the envelope for overlapping swallowtails). The outer 
envelope corresponds to # 7 = 90°, and moving inward the line type alternates as Q 1 decreases 
by 2°. 
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axis ratio q 



Fig. 6. — Similar to Figure 4, but for values of the shear angle near 6> 7 ~ 45°. Note that 
the vertical axes are log 7 instead of 7. For 9 y > 50.45° there are two closed envelopes. For 
50.45° > 9 y > 45° there is one closed envelope and one envelope that extends to 7 — > +00. 
For 45° > 9 1 > 39.55° there is only the envelope that extends to infinity. For # 7 < 39.55° 
there is no envelope. (Also see Figure 12 in Appendix A.) 




Fig. 7. — The dependence of swallowtails on the galaxy core radius s. All models have 
q = 0.5 and 7 = 0.22. The top panels show models with Q 1 = 88°, while the bottom panels 
show models with 9 1 = 90°. For comparison, Figures 2b and 2c show the same models in 
the limit s = 0. 
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Fig. 8. — The caustics for a galaxy with q = 0.5 perturbed by an external shear (top panels), 
and by a singular isothermal sphere at angle d 1 = 90° and distance d (bottom panels). The 
strength of the perturbation is given by 7. Each frame is 0.6 &sis,i on a side. In some cases 
with large 7 and small d, the caustics of the main lens galaxy and the perturber merge and 
become larger than the frames. 
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Fig. 9. — Similar to Figure 8, but for 6> 7 = 88°. Each frame is 0.6&sis,i on a side. 
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Fig. 10. — Differential cross sections (in units of 6| IS ) vs. logarithm of magnification A for 
2-, 4-, 6-, and 8-image lenses produced by galaxy+shear lens models. The cross-sections are 
computed for a point source in a square with side- length of 5 &sis- The lens galaxy axis ratio 
is q = 0.5, the shear amplitude is 7 = 0.22, and the shear direction is indicated; the cases 
correspond to those shown in Figure 2. 
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Fig. 11. — Sample images for a source with finite radius -R src . The lens model and source 
position are the same as in Figure 2b. The images are distorted tangentially relative to the 
lens galaxy, and as the source gets bigger they smear out into an Einstein ring. 
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axis ratio q 

Fig. 12. — The (q, # 7 ) plane for galaxy+shear models. The inset shows a closeup of the 
region with q small and # 7 ~ 45°. The solid curves show solutions to eq. (|A7|), while 
the dotted curves show solutions to eq. flA8f) . The solid curves intersect at 9 1 = 45° and 

q = ^7/(493 + 90^) = 0.0843. The dotted curves meet the solid curves at q = 8 - 3^ = 
0.0627 and c 7 = ±l/\/28 or 9 y = 39.55° and 50.45°. The curves divide the plane into 
nine regions. Each region has a distinct set of envelopes containing swallowtail models, as 
explained in Appendix A. 



